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Abstract 

Neutrino flavor oscillations are discussed in terms of two coupled Dirac fields. 
The interacting Lagrangian is diagonalized to obtain the exact eigenvalues and 
eigenf unctions. Flavor wave functions are then derived directly from the quan- 
tized neutrino fields. Probability densities obtained by squaring these wave 
functions upon taking into account the neutrino chirality are in agreement 
with the standard neutrino oscillation probabilities. 
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I. INTRODUCTION 



The theory and phenomenology of neutrino flavor mixing has been extensively studied 
mainly in the framework of quantum mechanics |l]-|3[]. Only very recently, a quantum field 
theoretical analysis of flavor mixing has been considered using the LSZ formalism 0]. 

Given the conflicting experimental results which have been obtained from a variety of 
neutrino observations, it is important to investigate the problem of neutrino propagation in 
the general context of field theory to reassure ourselves that there are indeed no differences 
between these results and standard ones obtained from the quantum mechanical treatment. 
Otherwise said, we need to firmly establish the approximations necessary to derive the 
expressions used phenomenologically. It is perhaps not superfluous to point out that in the 
Dirac theory, the contribution of negative energy states becomes substantial, an aspect of 
the problem totally absent in the quantum mechanical treatment. 

Our discussion in the present paper will be through the following Lagrangian 

£ = V> e («7 ■ d - m e )ip e + ^(ry ■ d - - 5(Vv0/* + Vy0e)> ( L1 ) 

consisting of two coupled Dirac neutrino fields ipe, with masses m e and m^. The inter- 
action is provided through a lepton number violating term with a coupling constant 5. The 
model allows for exact diagonalization. Neutrino and anti-neutrino flavor wave functions 
can be directly obtained as matrix elements of the quantized neutrino fields. 

The fully interacting Lagrangian should also include the weak interaction term. It is 
ignored here (as is done usually), since we are interested in free propagation and possible 
oscillation in flavor alone. For its production, we simply assume that the neutrino is created 
through some weak interaction process. Such a breakup is of course an approximation. 
The effect of parity non-conservation due to the weak interaction term on the other hand, is 
taken into account by considering the left-handed (right-handed) components of the neutrino 
(anti-neutrino) flavor wave functions. 

The model discussed here represents an example (albeit simple) of an interacting field 
theory which is exactly solvable. This may be of some help in elucidating the properties 
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of interacting theories, which can be normally studied only in perturbation theory. For 
example, the interaction between the electron and muon neutrino fields produces a change 
in their masses: the "experimental" masses mi, m<i depend on the "bare" muon and electron 
neutrino masses m e , m M and the coupling constant S. 

The paper is organized as follows. In Sec. 2, the standard neutrino oscillation phe- 
nomenology is reviewed. In Sec. 3, the two coupled Dirac equations, obtained from the 
Lagrangian (1.1) are solved. The electron and muon neutrino fields are then quantized in 
terms of the two free uncoupled fields, which diagonalize Eq. (1.1), as described in Sec. 4. 
The total conserved charge is the sum of the electron and muon flavor charges, which are 
not conserved separately. The last section closes with some concluding remarks. 

II. STANDARD TREATMENT 

We will critically review the standard quantum mechanical treatment of neutrino oscilla- 
tions to bring out the essential approximations made implicitly and the boundary conditions 
which are imposed. 

A state vector \ip > is introduced as a linear combination of the flavor eigenstates |e > 
and \ji > (assuming just two flavors) 

\ip >= C e \e > +C M |/i >, (2.1a) 

^>=(^)> (2.16) 

with C e =< e\ip >, =< fi\ip > and 

|C e | 2 + |C M | 2 = l. (2.1c) 

C e and then become the amplitudes for detecting an electron neutrino and a muon 
neutrino respectively. 

To derive the time evolution of the coefficients C e (t) and C M (i), the state vector \ip > is 
written as a superposition of the energy (mass) eigenstates \ui > and |z/ 2 > 

\iP >= C 1 \v 1 > +C2W2 >, (2.2a) 



' c 

\1>>=\ 1 I , (2-25) 



2 



(2.5) 



with Ci =< h>i\ip >, C 2 =< u 2 \ip > and 

|Ci| 2 + \C 2 \ 2 = 1, (2.2c) 

where Ci and C 2 are the amplitudes for finding the neutrino in the energy states E\ and E 2 
respectively. These coefficients evolve in time as 

d(f) = C 1 (0)e- iElt , C 2 (£) = C 2 (0)e-^ 2t . (2.3) 

Introducing the rotation matrix between flavor and mass eigenstates 

K>\ / cos(9) -sin{9)\ ( \e > \ 

(2-4) 

|i/ 2 > / \sin{9) cos{9) J \\fJ,>J 
it is easy to see that the following relation between the energy and flavor amplitudes holds 

'Ci(t)\ /cos(6) -sin(9)\ ( C e (t) 
C 2 (t)J \sin(9) cos(9) ) \C^(t) 
Hence, the time evolution of the coefficients C e and C M is given by 
C e {t)\ / cos(6) sin{9)\ /d(0)e- iElt ' 
C^t)) \-sin(9) cos(9)J \C 2 (0)e- iE ^ : 
In Eq. (2.6), the boundary condition that at production we have only a given flavor must 
be imposed. Suppose for example that at t = a muon neutrino is produced, i.e. 

C„(0) = 1, C c (0)=0. (2.7) 

/,From Eq. (2.5) at £ = we obtain 

Ci(0) = -sm(0), C 2 (0) = cos(0). (2.8) 

The time evolution of the flavor amplitudes is obtained by substituting Eq. (2.8) in Eq. 
(2.6) 

C e (t) = sin(9)cos(9)(e- iE2t - e^ 1 *), (2.9a) 



(2.6) 



C^t) = sin 2 (9)e~ iElt + cos 2 (6)e- im . (2.9b) 

Space and therefore momentum is introduced by assuming in Eqs. (2.9) E 2 = mj+p 2 , 
E\ = m 2 . + p 2 and x = t. The probability of finding a given flavor is obtained by squaring 
Eqs. (2.9) 

\C e \ 2 = sin 2 (26)sin 2 [(E 2 - E x )t/2] ~ sin 2 (26)sin 2 

\C^\ 2 = 1 - sin 2 (29)sin 2 [(E 2 - Ejt/2] ~ 1 - sin 2 (29)sin 2 

The assumption that the muon neutrino is created with a definite momentum p is only 
an approximation as has been pointed out previously |5|-f|. It is in contradiction with four 
momentum conservation, for example for the reaction tt — > fiu. Each of the possible energy 
eigenstates has a somewhat different momentum p^. In the rest frame of the pion, energy 
conservation dictates that (i= 1, 2) 

M n = ^jMl + p 2 + ^m 2 + V 2 . (2.11) 

Therefore, if we introduce momentum, i.e. space in Eqs. (2.9) we should write 

C e {x,t) = sin(6)cos(6)(e ipiX e- tElt - e ip2X e~ lE2t ) , (2.12a) 

C M (x,t) = sin 2 (6)e ipiX e- iElt + cos 2 (6)e ip2X e- iE2t . (2.126) 

In the relativistic approximation x ~ t the squared moduli of the amplitudes C e (x,t), 
C^(x,t) reduce to the standard neutrino oscillation probabilities given by Eqs. (2.10), as 
described in Ref. ||. The semiclassical approximation that a neutrino moves at a velocity 
close to c (assuming very small neutrino masses) "on a classical path" (x = ct) can be a 
good approximation if the neutrino travels over a macroscopic distance. 

Hence, if a muon neutrino is produced at the space-time (x ~ ct — 0), the probability 
of observing an electron neutrino is maximum for those space-time points x ~ ct at which 
\C e \ 2 = 1. If some matter is present at these points, processes such as v e + n — > p + e~, will 
occur, but not processes of the type v^ + n — > p + fi~ indicating therefore flavor oscillations. 
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(mf — m 2 )x 
2p 



(2.10a) 



(mi — m 2 )x 
2p 



(2.106) 



III. FIELD THEORETICAL DISCUSSION 



Since neutrinos are relativistic particles of spin 1/2, it is important to derive a relativistic 
equation of motion which can describe such flavor mixing. Energy eigenfunctions can be 
derived from this equation and one proper way to deal with states of negative energies is to 
quantize the field. As stated in the Introduction, we will consider the following interacting 
Lagrangian 

C = 4> e (i>y ■ d - m e )^ e + ^(ry • d - m M )Vv - <KVWv + ^M- (3-1) 

The parameter 5 is an extra mass (energy) related to the small amplitude that a neutrino 
can flip flavor. The following two coupled Dirac equations 

(i-f-d- m e )^ e - Si/jf, = 0, (3.2) 

(i-f-d- m^jVv - 5ip e = 0, (3.3) 

describe a neutrino which has some probability to flip flavor. If no flavor flips were possible 
(5 = 0), the rest energies (masses) of the system would be m e and m^, possibly equal to zero. 
However, since there is some amplitude that a neutrino, which is produced as an electron 
neutrino becomes later a muon neutrino, the possible rest energies of the system are not 
simply m e and m M , but are functions of the flipping energy. 
It is easy to see that the conserved current is 

^ = 4>vl^v = i>el^e + = J? + J£- (3-4) 

Thus, the separate electron and muon flavor currents are not conserved, only their sum is. 

In order to determine the energy eigenvalues and eigenfunctions of the system of equa- 
tions (3.2) and (3.3) we consider the ansatz 

^ = ae- iPx , (3.5) 

= be~ tPx , (3.6) 
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where P is the four-momentum P = (E,p), which is unknown and is to be determined so 
that the system of differential equations (3.5) and (3.6) is satisfied. The coefficients a and b 
are Dirac spinors, which can be written as 

/Xi\ 

a = , (3.7a) 

\X2/ 

6= (r)' {s ' 7b) 

where Xi,2 and y?i,2 are two component vectors. Substituting Eqs. (3.5), (3.6) into Eqs. 
(3.2), (3.3), we obtain the system of linear homogeneous equations 

Exi = (^-PX2 + m e xi + S(p 1 , Ex2 = & ■ PXi ~ ™-eX2 - 5(p 2 , (3.8a) 

E(fi = a ■ py?2 + m fl tp 1 + Sxi, Eip 2 = a ■ p</?i - ~ $X2, (3.86) 

where a are the Pauli matrices. 

The system of Eqs. (3.8) admits non trivial solutions only if 

E A - E 2 (2p 2 + 25 2 + m 2 e + ml) + p 4 + 5 4 + 

p 2 (25 2 + m 2 e + ml) + n^m 2 - 25 2 m e m^ = 0. (3.9) 
Solving Eq. (3.9), we obtain (p = |p|) 

E 1>2 = ±yjp> + ml 2 , (3.10) 



with m 12 given by 



and 



m 1>2 = -[(m e + m M ) ± R], (3.11) 



R= ^ mii - me ) 2 + A5 2 . (3.12) 

Therefore, while in the free Dirac equation there are two energies (one positive and one 
negative) for every possible value of the momentum p, for a system of two coupled Dirac 
equations, for every possible value of the momentum, there are four possible values of the 
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energy, two positive and two negative, due to the possibility of (flavor) oscillations. Also, 
because there is some chance that the neutrino can flip flavor, the rest energies of the electron 
and muon neutrino system are not simply m e , but are given by Eq. (3.11). 
The solutions of Eqs. (3.8) can be written as 

S(a ■ p)(ro e + m^) _ 5{E 2 - p 2 - 5 2 - E(m e + r?y) + m e m M ) 

Xi - q V?2, X2- Q ^2, 

(a-p)(p 2 + m 2 e + 5 2 -E 2 ) 

<Pi = q V2, (3.13) 

with 

Q = —E 3 + E 2 m^ + E(5 2 +p 2 + m 2 e ) + 5 2 m e - p 2 m^ - m 2 e m^. (3.14) 

As stated above, for a given value of the momentum p, there are four different energies 
±Ei t 2 and for each energy, there are two eigenfunctions with different (up and down) spins. 



Corresponding to the positive energy solution E 1 = \jm\ + p 2 , we have the following two 
solutions 

M*,t) = ^-^=M^P)e ip - x e-^\ (3.15) 

where s — 1, 2 is the spin index and 0i(s, p) is given by 

1 / «i(s,p) \ 
M*,P)= , h (3-16) 



with 



and Ui(s, p) is the Dirac spinor 



m„ — m P + R , 

M 1 = * e . (3.17) 

2o 



tti(s,p) = y/Ea +mi | |, (3.18a) 

Ei+mi A 

satisfying the Dirac equation (7^^ — mi)ui(s, p) = 0, and y^" 1 satisfies the normalization 
condition 

X tw X W = 1- (3-186) 



7 



For the other positive energy solution E 2 = \Jp 2 + mf, we have 



with 2 (s, p) given by 



and M 2 defined as 



M*,t) = ^-^=Ms,p)e^e-^\ (3.19) 



1 



M 2 = ^ . (3.21) 



/ «2(s,p) 

fc( a ,p) = -== , (3.20) 

/ 1 + M| \M 2 u 2 S)P 



We notice here that because M ± M 2 = —1 we can write 2 (s, p) in terms of Mi as 

1 



Mitt 2 (s,p) 

Us,p)= r== I |, (3-22) 

^l + Mf V -u 2 (s,p) 



and w 2 (s, p) is the Dirac spinor 



M 2 (S 



x w 

,p) = V^ + ^ ,J, (3-23) 



satisfying the Dirac equation (^p^ — m 2 )« 2 (s, p) = 0. 

Similarly for the solutions of negative energies — E i2 we have the following eigenfunctions: 



M*> t) = ^^&(*, p^V**', (3.24) 



with 03 (s, p) given by 



1 /-M 2 ui(s,p)\ ! / vi(s,p) 

3 ( S) p) = = . (3.25) 

Vl + MlV Vl ( s ,p) / yl + M? \M 1?)l ( S) p) 



and 



g P v (a) 
A 



Wl ( S ,p) = v^iT^r i+ ";\ , (3.26) 



with 04 (s, p) given by 



^ 4 (x, f) = -L-J^ 4 ( a , p)e~^e + ^\ (3.27) 



x / Miv 2 (s,p) 
Ms,p)= I , (3.28) 

VI + Ml \ -v 2 (s,p) 



and 



t^( S , p) = 4E 2 + mi f ^ J . (3.29) 



IV. FIELD QUANTIZATION, ANTI-COMMUTATION RELATIONS, AND WAVE 

FUNCTIONS 



We introduce the matrix U 



Mi 



jj = I V 1+M ? y/ 1+M ? 

M± 1 

to make the transformation on the field 

f^e\ (<I>1 
WJ V02, 

which uncouples the interacting Lagrangian given by Eq. (3.1) 



(4.1) 



(4.2) 



£d = 0i(«7 ' ^ - ™i)0i + 02^7 • d - m 2 )4>2- 



(4.3) 



Following relations are useful to see that U uncouples Eq. (3.1) 



6 (Mi + M 2 ) = m M - m e ; R = S(M 1 - M 2 ); MiM 2 = -1. 



(4.4) 



Therefore the fields 0i and 2 describe the "normal modes". To quantize ip u , we expand 
ip e and ^ M in terms of the normal modes (energy eigenfunctions) found in Sec. 3 

1 1 



1p e (x,t) 



EE 



bi(s, p)™e^ + M M s, p)WW e -« 



-tpx 



v^^l + Mf p 



EE 



v 7 ^ 



-e-^-6 2 ( S ,p)- 



M I 4( S ,p)^--4( S ,P)^e-« 



tpx 



e* p ' x + 



(4.5a) 



e 4p ' x + 



(4.56) 



where the number operators 6j and dj (i=l,2) satisfy the anti-commutation relations 



{&i(s,p), &t(s',p')} = SijSrfSss,, 
{di(s,p), d](s',p')} = SijSpp'Sss', 



(4.6a) 
(4.66) 



and all the other anti-commutators are zero. 61,2(5, p) annihilates the normal mode of 
positive energy E 1>2 and spin s; d[ 2 creates the anti-normal mode of positive energy E 12 
and spin s. The vacuum state is defined by 

6 i |0>=d i |0>=0. (4.7) 

The total charge operator is 

Q = Qe + Q„ = J ^vd 3 X = J (^Ve + ^»)d 3 * = 

EE[ 6 i( s 'P) 6 i( s 'P) P)M«>P) -4(s,P)di(s,P) -4(s,p)4(s,p)]. (4.8) 

p * 

The following relations hold 

[Q, 6}(s, p)] = 6j(s, p) [Q, d\(s, p)] = p). (4.9) 

Hence, for a given value of the momentum p and spin s, there are four possible normal mode 
states 

b[(s,p)\0>=\l pa >, 6 2 (s,p)|0>= |2 PS >, (4.10a) 
4(*,P)|0 >= I - l ps >, 4(s,P)|0 >= I - 2 ps > . (4.106) 
These states differ for the charge ( ±1, ±2) 

Q|l ps >= |l ps >, g|2 ps >= |2 ps >, (4.11a) 

Q\ - l ps >= -| - l ps >, Q\ -2 ps >= -I -2 ps >, (4.116) 

The above states allow us to construct wave functions in space-time. For example, the 
wave function associated with the state |l ps > is 



Vv(x,t) 



> C (X,*)\ /<0|4(x,t)|l ps > 

,W(x,o/ \<o|^(x,o|i ps >. 



52 J ^«'(^)^'' < 4 - 12 > 
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This being a plane wave, gives a stationary probability of finding a neutrino at a given 
space-time point. However in any location inside the volume V there is a probability equal 
to ( i+ M i ) of finding the neutrino in the electron flavor and probability equal to ( ) of 
finding it in the muon flavor. Similar considerations can be applied for the other three states 
|2 ps >, | — l ps >, and | — 2 ps >. Therefore, these states represent states of mixed flavor 
at any given space-time point. To obtain states for which we have only one flavor at a 
given space-time point we need a superposition of states. A general state of positive charge, 
momentum p and spin s is given by 

|0 + >= [Ab\(s,p) + Bbl(s,p)}\0 >, (4.13) 

where A and B specify the amount of each normal mode state of positive energy present in 
the state \<f>+ >■ 

Similarly, a general state of negative charge, momentum p and spin s is given by 

|0_ >= [Cd{(s, p) +Dd\(s, p)]|0 >, (4.14) 

where C and D specify the amount of each normal mode of negative energy present in the 
state |0- >. The matrix element 

1 1 



< O|^ e (x,t)|0+ >= V> e (x,f) = 



W J I + Ml 



(4.15a) 

gives the probability amplitude of finding a neutrino of momentum p and spin s at the 
space-time point (x, t) with the electron flavor. In the same way, the matrix element 



< o|Vv(x,t)|0 + >=Vv( x >0 = 1 ' 



MiA M^9) e - iElt _ B Ms^ e _ lE2t 



(4.156) 

is the probability amplitude for the muon flavor. 

Similarly the matrix elements < O|-0 e (x, > and < O|'0 (U (x, t)\<f>- > give the proba- 
bility amplitudes for finding an electron anti-neutrino flavor and a muon anti-neutrino flavor 
respectively. 
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To take into account the fact that neutrinos (anti-neutrinos) are created with negative 
(positive) chiralities, we define the "observable wave functions" as 

^ eL (x,t) = (l- 75 )^e(x,t), VVL(x,t) = (l- 75 )^(x,f) (4.16) 

where i/) e (x., t) and ^(x, t) are given by Eqs. (4.15). 

If p is for example along the z-axis, we can choose x^ as 

X (1) =n, X i2) =(°), (4.17) 



and the left-handed Dirac spinor is 
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y(2) \ / y(2) 



„ lfas = 2 ) = (l- 75 )V¥+^ =V¥+^(1 + ^-) . (4.18) 

V E+mX. / \ X / 

In the limit E ^> m, (1 — 75 )-u(p, s = 1) ~ 0. 
For p along an arbitrary direction n 

P = p[e x sin{6)cos((f)) + e y sin(0)sin(<f>) + e z cos(9)}, (4.19) 

we can choose the spinors x^ as 

/e-i*cos(§)\ /-e-^sm(§)\ 

* (1)= ' * (2)= > ■ (4 - 20) 

V e5*sin(|) / V e2*cos(f) / 

and the left-handed spinor is still given by Eq. (4.18). 

Hence, the observable flavor neutrino wave functions are 

,t, , (x t ) - 1 e ""~ \A VEl + mi a + p )c- E " 



(2) 



12 



The coefficients A and B in Eqs. (4.21) are determined through the initial boundary 
conditions. Suppose that at t — 



^ L (x,t = 0) = 0, 



(4.22) 



we have only the electron flavor present. The other one is obtained by the normalization 
condition 

/ rf 3 x|^ eL (x,t = 0)| 2 = 1. (4.23) 
J v 

By imposing the boundary conditions given by Eq. (4.22) and Eq. (4.23) we obtain the 
following flavor wave functions 



^eL(x,t) 



3 !P'X 



V \ + M%y/V 



-iEit 



+ M(e 



2iE 2 U 



X 



(2) 



ipx jur 1 

\c~ iElt - e iE2t 



X 



X 

(2) 



(2) 



(4.24a) 



(4.246) 



Hence the probability densities of finding the electron and muon neutrino flavor are given 
respectively by 



Pe(t) = 



V 



1 - 



2Mi 
1 + Mj 2 



(£ 2 - Ejt 



1 / 2Mi 



sin' 



(E 2 - 



(4.25a) 



(4.256) 



V\1 + M?J 2 
The coefficient [2Mi/(l + M 2 )] 2 is equivalent to sin 2 (29) in Eqs. (2.10). Therefore the 
field theory treatment reduces to the standard quantum mechanical treatment described in 
Sec. 2. 

As another example, we consider the case in which the neutrino eigenfunctions of different 
masses have different momenta pi and p 2 with 



Ei = Vpi + ^i, E 2 = ypi + 



(4.26) 



The neutrino flavor wave functions Vvl(x, t) are therefore 

in 



■(x,t) = * M ! f9 4[e"' El V pi ' x - e lB2 V P2 ' x ] f 



(2) 



-x 



(2) 



(4.27a) 
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111 / x {2) \ 

where we have assumed that at a given space-time point (x = 0, t = 0), we have only the 
muon flavor. 

The probability densities of finding at the space-time point (x, t) the electron and muon 
neutrino flavors reduce to Eqs. (2.10) in the relativistic approximation |x| ~ t. 

V. CONCLUSIONS 

We have discussed an explicit model of neutrino flavor mixing in the framework of quan- 
tum field theory. In this model, the equations of motion for the interacting fields are solved 
directly and the system is diagonalized in terms of the two uncoupled free fields <f>\ and 02 
of mass mi and rri2 respectively. We notice here that because we can directly diagonalize 
the Lagrangian we do not need to write the interacting fields in terms of the free asymptotic 
fields ipoe, V'o/i °f m ass m e and respectively. We have also derived neutrino flavor wave 
functions in such a way that the total flavor charge is constant. The probability densities, 
derived from these wave functions, are in agreement with the standard neutrino oscillation 
probabilities, if we take into account the neutrino chirality. 

Also, since explicit plane wave solutions for all normal modes have been obtained, wave 
packets corresponding to these can be constructed via standard techniques ||. 
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